Coherent trajectory through the normal-to-superconducting transition reveals 
ultrafast vortex dynamics in a superconductor. 
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We report on measurements of the trajectory of the superconducting order parameter ip(t) through 
the normal-to-superconductor transition in Lai.gSro.iCuCU under ultrafast quench conditions up to 
4 x 10 14 K/s. Using a 3-pulse technique, we isolate the trajectory of ip(t) from single-particle pro- 
cesses on short timescales. Detailed modeling of the system trajectory using hot electron energy 
relaxation, time-dependent Ginzburg-Landau theory and vortex creation according to the Kibble- 
Zurek mechanism describes the experimental data in quantitative detail, significantly advancing our 
understanding of the behaviour of the normal-to-superconductor transition under highly nonequi- 
librium conditions. The recovery of takes place on a timescale of 20 to 60 ps, depending on the 
excitation energy density, while the vortex dynamics is apparent on a timescale of 10 ~ 30 ps and 
is visible only for the fastest quench rates. 



What happens when a system goes through a 
symmetry-breaking transition depends crucially on how 
fast the transition occurs [TJ [2]. If the quench process is 
sufficiently fast large numbers of topological defects are 
created and the nature of the final state can be substan- 
tially different than if the transition takes place slowly. 
The study of the time evolution of complex systems 
through symmetry breaking transitions (SBT) is thus of 
great fundamental interest in different areas of physics 0- 
4J and beyond. The diversity of microscopic vacua in 
condensed matter systems makes them excellent candi- 
dates for the experimental study of new phenomena[5j, 
also posing new challenges for the general theory of 
SBTs. Of particular general interest is the normal- 
to-superconducting (N— >*S) state transition in which a 
Lorentz non- invariant system breaks gauge invar iance[6j. 
Under near-equilibrium conditions the N^S transition 
occurs when the thermal fluctuation timescale rp be- 
comes comparable to the "Ginzburg-Landau time" which 
for typical high temperature superconductors is tql = 
7r/h(T - T c )~ 10- 12 - 10- 13 s[2j |5j E]. The intrinsic 
collective system relaxation time tql defines the relevant 
timescale and experimental resolution necessary to in- 
vestigate system trajectories through the transition. Till 
now, the shortest quench times through the transition 
were r q ~ 10 -8 s[9j, which were relatively long compared 
to tql , so observations were focused on statistical counts 
of vortices in the aftermath of the transition as a func- 
tion of the quench rate 7 q = l/r q . Femtosecond laser 
spectroscopy offers sufficient time resolution for ultra- 
fast quench studies and pump-probe transient reflectiv- 
ity experiments have recently elucidated the destruction 
of superconductivity (i.e. the S— ?>N transition) [TDHH] . 
However, standard pump-probe techniques cannot dis- 
tinguish the trajectory of the order parameter (OP) ip(t) 
from photo-excited hot carrier relaxation, quasiparticle 



(QP) recombination dynamics across the superconduct- 
ing gap and the relaxation of pseudogap (PG) excita- 
tions, all of which appear simultaneously in the time- 
resolved response [T3HT8j . 

To solve this problem we devised a 3 pulse tech- 
nique shown schematically in Fig. 1 which allows us 
to measure in real time the superconducting order pa- 
rameter ip(t) through the transition in a Lai.gSr .iCuO4 
(LSCO) single crystal with ultrahigh quench rates of up 
to 7 q = l/r q = (dT/d^T' 1 ~ 4 x 10 14 s" 1 . Further- 
more, we can vary j q by changing the laser excitation 
energy. We compare the data with predictions of the 
Kibble-Zurek (K-Z) mechanism of vortex creation, find- 
ing evidence for ultrafast vortex dynamics in the rapid 
quench regime when r q ~ tgl[1> 20J. Lai.gSro.iCu04 is 
chosen because it is one of the best studied superconduc- 
tors particularly by femto-optical methods. [TO} [T4| [27] 

With our technique the superconducting state is first 
rapidly nonthermaly[10j destroyed in ~1 ps by an intense 
"destruction" (D) laser pulse creating a pancake-like nor- 
mal region defined by the laser beam profile with the 
FWHMdiameter d = 60 /im and the optical penetration 
length A p = 150 nm (at 800nm). We then monitor the 
recovery of the OP by means of the transient reflectiv- 
ity AR(At2s)/R excited with a weaker pump (P) pulse 
and probed with even weaker probe (p) pulse as a func- 
tion of time delay Ati_2 between the D and P pulse (see 
Fig. 1).[H] The D pulse fluence is adjusted to achieve 
different quench rates. 

After destruction, an OP front propagates from 
the bulk towards the surface on a timescale ~ 
A p /^ ~ A p tgl/£± ~ 10 3 tql, where is the OP front 
velocity [32J and £_l is the out-of-plane SC correlation 
length. T^jj is much longer than both, the quench time 
r q and the Zurek time rz ~ ^Jt q T gl . As a consequence 
the normal region is supercooled and we can expect vor- 
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Figure 1: a) A schematic diagram of the pulse sequence. The 
time delays Ati_2, Ati-3 and At2-3 refer to delays between 
the D pulse (blue) and the P pulse (red), the D pulse and the 
probe pulse (green) and between the P pulse and the probe 
pulse, respectively. The temperature front where T e = T c 
and the S/N phase boundary move with velocities vt and 
respectively. Vortices are created in the wake of the temper- 
ature front, b) The system trajectory (depicted by the silver 
ball) in a temporally evolving potential. In the rapid quench 
scenario (A), the potential changes faster than the system can 
follow. The opposite is true in the slow quench scenario (B). 



tex formation according to the K-Z scenario. [TJ 15] 19} [201 
[22] Depending on the quench conditions, the propagat- 
ing velocity of the temperature T — T c front can 
be either faster or slower than the critical velocity [20J 
v cr i t « yj tgl / 7q ^ 10 5 cm/s. In the ultrafast quench 
limit when vt > ^crit? we expect vortex formation, while 
in the slow quench limit (vt < v cr n), vortex formation 
becomes suppressed [20) [22]. 

The time-evolution of the electronic and lattice tem- 
peratures T e (t) and T^(t) respectively, which govern the 
quench process, are calculated using the standard two- 
temperature model[231 leT e ^f = -Jh(T e - T L ) + 
E(t) and C L (T)^ = 7L (T e — Tl), where E(t) is a func- 
tion (Gaussian) describing the energy density per unit 
time supplied by the laser pulse [33J. For maximum 
accuracy we used the temperature-dependent thermal 
constants 7 e (T) and Cl(T) from experimental data[25] 
and the measured energy loss rate 7l — 340 K/ps for 
Lai. 9 Sr .iCuO 4 [24J. The dependence of T e (t) and T L (t) 
on the D pulse fluence, J 7 ^, is shown in Fig. 2a) for the 
first 2.5 picoseconds. The red contour lines show when 
either T e or Tl crosses T c , which is the critical point of 
the quench. 

For low J~y> < 15/iJ/cm 2 , we have a rapid quench sce- 
nario as T e rapidly crosses T c within a few hundred fem- 
toseconds, the quench rate being governed by the elec- 
tronic energy relaxation rate 7 q c± j e ~ 4xl0 14 K/s. For 
J~d >15 /iJ/cm 2 however, the lattice is heated above T c 
by the D pulse, so the cooling rate through T c is gov- 
erned by diffusion processes on timescales longer than 
~ 30 ps. Thus by changing Fy> we can effectively vary 
the regime from rapid to slow quench, over a very wide 
range of quench times from r q = 100 fs to hundreds of 
picoseconds. 




Figure 2: The time evolution of the electronic and lattice 
temperatures T e and Tl respectively. The ball signifies the 
trajectory of T e (t) for a D pulse fluence Jr> — 38 /xJ/cm 2 . 
The blue lines are the trajectories of T e (t) corresponding to 
the D pulse fluences used in the experiments. The red lines 
signify the time when T e and Tl cross the critical temperature 
T c = 30 K. Rapid quench occurs for Tt> < 15 /xJ/cm 2 . 



The transient reflectivity amplitude in the supercon- 
ducting state is proportional to the square of the order 
parameter A s = AR/R oc 1^(^1-2) | 2 [E] 5 so to detect the 
trajectory ip(t) through T c , we measure the transient re- 
flectivity by the Pump-probe technique at different time 
delays after the D pulse. The experimental data on the 
transient reflectivity AR(At2s)/R at different delays 
Ati_2 during the system recovery is shown in Fig. 3a). 
Two distinct and easily identifiable contributions are ob- 
served: one has a peak around 0.2 ps which has been 
previously identified to be due to pseudogap excitations, 
and the other is the QP recombination signal of the su- 
perconducting condensate with a peak near 2 ps, extend- 
ing to tens of psp~4j[27]. We observe a clear suppression 
of the QP peak at short delay Ati-2 and a gradual emer- 
gence of the SC state with increasing delay. A distinct 
feature in Fig. 3a) is the convergence of all QP curves 
after a delay At2s = ^th — 25 ps which corresponds to 
the complete thermalization of all degrees of freedom. 

For further analysis, we subtract the PG response from 
the data[34j and plot A s as a function Ati_2, thus obtain- 
ing the trajectory of \ip(t)\ 2 shown by the square black 
dots in Fig. 3b). From the exponential fits of the decay of 
AR/R(t2s) we also obtain the QP relaxation time tqp 
as a function of At 1-2 which is shown in Fig 3b), plotted 
as a relaxation rate 1 / tqp . We observe that 1 / tqp shows 
a similar trajectory as the order parameter; indeed su- 
perimposing |^(ti_2)| on the graph of 1/tqp(£i_2) (solid 
line in Fig. 3b)), we see that the QP recombination 
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Figure 3: The transient reflectivity data AR/R for 
Lai.gSro.iCuCU at 4 K as a function of At2-3, for different 
At i-2. The D pulse fluence is Tt> — 12/i J/ cm 2 , which is 
approximately three times above the destruction threshold 
(J-d — 4/i J /cm 2 ) [TO]. The red line indicates the pseudogap 
signal measured above T c . The recovery time of the SC state 
Tree is indicated, b) The same data as in a) after subtrac- 
tion of the PG signal plotted as a function of At 1-3 . The 
black squares are the fitted maximum amplitudes of A s and 
give ^(t)] 2 for each At 1-2. The blue squares are the QP re- 
laxation rate superimposed on the measured order parameter 
trajectory \ip(t)\ oc y/A^ as a function of delay Ati-2. 



rate closely follows the order parameter as a function of 
Ati_2, consistent with the fact that 1/tqp(£) is propor- 
tional to the superconducting gap A s (t), which is in turn 
proportional to \ip(t)\ [13J. 

Tuning Tr> significantly changes the trajectory of the 
system through the transition, and the SC state recovers 
much more slowly with increasing J 7 ^. Fig. 4a) shows 
the trajectory of |^| 2 for a number of different ^D-i^J 
For the lowest two values of J 7 ^, which correspond to the 
fastest quench, we see that destruction of the condensate 
is incomplete throughout the entire excitation volume, 
since |^| 2 is non-zero in the first 2 ps after the D pulse. 
With increasing J 7 ^ the recovery becomes slower. The 
initial recovery rate 1/t[ can be obtained from the ini- 
tial slope of |^(Ati_2)| 2 for small Ati_2 and is shown 
in Fig. 4b). Plotting the recovery time versus J 7 ^ in 
Fig 4c), we find a phenomenological r\ which increases 
approximately linearly with J 7 ^ before eventually satu- 
rating at high J~i> For comparison, we see that the QP 
recombination time tqp (shown on the same plot) first 
dramatically decreases below the destruction threshold 
fluence for J 7 ^ < 4/iJ/cm 2 , and then slightly increases 
with increasing J 7 ^ . This impressively demonstrates how 
well we can separate the order parameter recovery from 
tqp in these experiments, and more importantly, shows 
that the QP relaxation has very different systematics and 
plays no significant role in determining the recovery dy- 
namics. 

To analyze the trajectory of the system through the 
transition as accurately as possible, we use the full set 
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Figure 4: a) The trajectories of \i/j\ 2 as a function of At 1-2 
for different D pulse fluences Fy>. The lines are calculated 
trajectories from the TDGL theory. The insert shows the 
predicted ^(t)] 2 with (black) and without (orange) vortices 
at 7/jL J / cm 2 . The vertical scales are normalized to the equilib- 
rium I -00 1 2 in the absence of the D pulse, b) The relaxation of 
I -0| 2 in the first 5 ps for different fluences. c) The dependence 
of the recovery time T rec on pulse fluence (red dots). For 
comparison, the QP recombination time tqp is also shown 
for different Tr> (diamonds). The dashed line is the predicted 
tqp from the Rot hwarf- Taylor model[28, 29 which is valid 
only below the destruction threshold Ty> < 4/xJ/cm 2 . 



of TDGL equations where in the absence of an external 
field, the dimensionless order parameter ip(t) is described 
by [30]: 



~dt 



-a r (t,^-V#| 2 + V 2 V> + 77, (1) 



(2) 



where temporal and spatial coordinates are measured in 
units of tql and correlation length £, respectively. The 
only parameter in the equation is the ratio u = I ^ be- 
tween the longitudinal relaxation time and the transverse 
relaxation time. For our LSCO sample, tq « 1 ps, u « 5 
[8 J and £ « 2 nm. To initiate vortex formation arising 
from microscopic fluctuations, the customary Langevin 
noise term 77 is added |31| . We also take into account the 
inhomogeneous nature of light penetration into the sam- 
ple and thermal diffusion processes, which take place on 
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Figure 5: a) The predicted vortex structure of at dif- 
ferent delay times corresponding to rapid quench conditions 
(Jt> = 12 pJ/cm 2 ) showing the \ip(r, t)\ 2 in the time interval 
Ati_2 = - 25 ps in an area of 200 x 200 £ 2 . b) The dif- 
ference between the experimental trajectory and the fits cal- 
culated without Langevin noise: A|^| 2 = iV'Inovortex — l^lexp 
shows a clear suppression of |^| 2 around 10-20 ps for 5 < 
Jr> < 12 J/cm 2 . c) the predicted suppression of A \ip\ theory — 
l^lnovortex — 1^1 Cortex due to vortex formation around 10-20 ps 
introduced by fluctuations. 



timescales beyond 10-30 ps by including a solution to the 
heat diffusion equation 9Tdi g/ t?s) = D(T)V 2 T diff , where 
D(T) = k/pCi, is the diffusion coefficient, p is the den- 
sity, k is the thermal conductivity. Solving for diffusion 
into the sample, we obtain a quench function applicable 
for long times. A simplified quench function can now be 
introduced in the TDGL equations: 



a T (t,z) 



-1 + ?7L 



4Dt 



2 /K P -t/n 



(3) 



where the second term is a simplified solution to the dif- 
fusion equation in the limit of negligible lateral diffusion, 
and the last term approximates the trajectory of T e (t) 
shown in Fig. 2a). The fitting parameters are the two 
constants 7?l and r] e . In Fig. 4a) we first compare the 
fits of the calculated trajectories with the data using Eq. 
[2] with r] = 0. The fits are seen to be impressively good 
over the entire trajectory for high values of .Fd, but fail 
in a region around 10-20 ps for low J 7 ^. (Note that this 
discrepancy is extremely reproducible and well above the 
noise level) . It is not possible to remove this discrepancy 
by changing any of the parameters of the model fit. Con- 
sidering the low scatter of the data and the quality of the 
fits for higher J 7 ^ it is clear that the persistent suppres- 
sion of |?/;| 2 around 10-20 ps in the fast quench regime is 
a real effect. 

Next, we include n in the calculation, initiating vortex 
formation. The resulting vortex structure is shown in 
Fig. 5a) in a sequence of snapshots using the parameters 
for the trajectory at = 12/iJ/cm 2 . The suppression of 



\il)\ 2 by vortices calculated for different is compared 
with the experimentally observed effect in Fig. 5 b) and 
c) showing a remarkable agreement with the data around 
10-20 ps and J 7 ^ < 14/iJ/cm 2 Intrinsic vortex-antivortex 
creation and annihilation dynamics is predicted precisely 
in the region in time and pulse fluence J 7 ^ (i-e. quench 
rate) where the data show deviations from the homoge- 
neous trajectory. In contrast to previous studies [9j [21], 
the experimental vortex density (Fig. 5b) appears to be 
in reasonable quantitative agreement with the calculated 
K-Z vortex density ( Fig. 5c). 

To conclude, we have found that macroscopic phe- 
nomenological modeling can accurately describe the data 
very well without resort to microscopic theory provided 
we include vortex dynamics within the Kibble- Zurek sce- 
nario. The microscopic phenomenology may be thought 
to be included in the quench function, which includes the 
effect of the quasiparticle relaxation bottleneck appearing 
as a result of the formation of a gap in the superconduct- 
ing state. One may try and introduce modifications to 
the quench function, but these would have no transparent 
physical meaning within the present context. An alter- 
native microscopic description would require a complete 
microscopic theory. The fact that no external adjustable 
parameters - other than those of the noise - are needed 
to achieve this agreement is attributed to the great care 
taken to accurately model the quench process and inho- 
mogeneity. The vortex-antivortex pair generation rate 
according to the Kibble-Zurek scenario[20 ) 22J is propor- 
tional to T~ a , where a = 0.25 ~ 0.5(50], consistent with 
the model fits and the observed suppression of xp under 
ultrafast quench conditions. We have demonstrated that 
our new experimental technique opens the way to real- 
time studies of the order parameter trajectory and vortex 
dynamics in superconductors, giving information which 
is not available from other time resolved methods, and 
shows that modeling by use of TDGL equations can give 
a faithful description of the basic phenomena presented 
here for the first time. We end by noting that the signifi- 
cance of the dependence of t\ on J 7 ^ needs to be clarified, 
which can only be done by further systematic studies of 
other systems. 

We wish to acknowledge the useful discussion with 
T.W. Kibble regarding the importance of a variable 
quench rate experiment. 
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